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We carefully analyze the renormalization group equations in the type I + II seesaw scenario in the 
extended standard model (SM) and minimal supersymmetric standard model (MSSM). Furthermore, 
we present analytic formulae of the mixing angles and phases and discuss the RG effect on the 
different mixing parameters in the type II seesaw scenario. The renormalization group equations of 
the angles are independent of the Majorana phases. In addition, there is a contribution which is 
proportional to the mass squared difference for a hierarchical spectrum. This is in contrast to the 
inverse proportionality to the mass squared difference in the effective field theory case. 
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I. INTRODUCTION 

In the last years, the precision of the leptonic mix- 
ing parameters has been further increased [ll due to a 
number of different experiments [1 H 0, HH ■ The 
precision of measurements of leptonic mixing parameters 
will even increase in upcoming experiments. In the next- 
generation experiments, the mixing parameters will be 
measured on a 10 % level @. Thus we are entering an 
era of precision experiments in neutrino physics. 

Contrarily, there have been a lot of attempts to ex- 
plain the structure of the neutrino mass matrix {e.g. 
H [H EH)- However, most models use heavy particles 
which generate light neutrino masses effectively after de- 
coupling. Examples are right-handed neutrinos in the 
standard (type I) seesaw mechanism [IH, [TJ, [TJ, [H 



m v = -(m^T M 



(1) 



or a Higgs triplet in the type II seesaw mechanism [161 

MM 



vaYa 



(2) 



where «a is the vacuum expectation value (vev) of the 
Higgs triplet and Ya is the Yukawa coupling matrix of the 
vertex IA£. Clearly, the seesaw operates at high energy 
scales while its implications are measured by experiments 
at low scales. Therefore, the neutrino masses given by 
Eqs. (|1I2[) are subject to quantum corrections, i.e. they 
are modified by renormalization group (RG) running. 

The running of neutrino masses and leptonic mixing 
angles has been studied intensely in the literature. RG 
effects can be very large for a quasi-degenerate neutrino 
mass hierarchy and they can have interesting implica- 
tions for model building. For instance, leptonic mix- 
ing angles can be magnified [H El EH, [H El, the 
small mass splittings can be generated from exactly de- 
generate light neutrinos [H El El ES El El or bi- 
maximal mixing at high energy can be made compati- 
ble with low-energy experiments [H [HI, HI- On the 
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other hand, even rather small RG corrections are impor- 
tant in view of the precision era we are entering. For 
example, RG effects induce deviations from #13 = or 
maximal mixing 023 = 7r/4 [H, EH, [35| also for a hier- 
archical spectrum, as well as from other symmetries, like 
quark-lepton-complementa rity l36l . [37t l38l El , tribimax- 
imal mixing 3jl 4(| [4l|, [42|, [43|, |44| | and other special con- 
figurations J49,j4fj. Threshold corrections can yield larg e 
effects [II [MiM El EI El El 51 B 

These studies have been done in the effective th eory 
of Majorana neutrino masses [H E3, El, [H [H, l55L l5al . 
There are also studies in the standard seesaw case [52|, [57| 
and in the case of Dirac neutrinos [HEl which also show 
significant RG effects which can become comparable to 
the precision of experimental data. Therefore the RG 
effects have to be considered in model building in order 
to be able to compare predictions to experimental data. 

In this paper we present the RG equations in the t ype 
II seesaw scenario [75[ in the standard model (SM) [6fJ 
and minimal supersymmetric standard model (MSSM). 
We derive analytic formulae which allow to understand 
the running of the neutrino parameters above the thresh- 
old of the Higgs triplet. Furthermore, we extend the 
software package REAP/MixingParameter Tools [7(| by 
an Higgs triplet for analyzing the RG evolution numeri- 
cally. A similar calculation has been done by Chao and 
Zhang on the renormalization of the SM extended by a 
Higgs triplet [rJil ]. The two calculations differ in several 
terms (77j |. 

The paper is organized as follows: In Sec. [Hi we present 
the Lagrangian of the type II seesaw model and give the 
tree-level matching conditions. Furthermore in Sec. IIIII 
we show all new wave function renormalization factors 
and counterterms and point out the differences to the 
work of Chao and Zhang [6l|. The RG equations are 
shown in Sec. IIVI In Sec. [V] the additional terms in 
the superpotential and the RG equations in the MSSM 
are presented [zl]. Sec. [VI] is dedicated to the analytic 
understanding of RG effects in the type II seesaw case 
(only a Higgs triplet) and Sec. IVIII gives a glimpse on the 
full type I + II seesaw case. Finally, Sec. IVIIII contains 
our conclusions. 
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II. TYPE II SEESAW LAGRANGIAN 

In the following, wc consider the SM extended by right- 
handed neutrinos ~ (1, 0)su(2)xU(i) an< ^ a charged 
Higgs triplet A - (3, l) SU (2)xU(i)> 



A+/V2 A++ \ 
A -A+/72 ) 



The Lagrangian is given by 



where the individual parts are defined by 



(3) 



(4) 



i?A =tr [(£> M A) f £>"a] - r(A, <t>) 

-^=(Y A ) fg e T L f aa 2 Aei + h.c. 

f(A, 0) =M| tr (At A) + ^ (tr A* A) 2 

+ ^ [(tr At A) 2 - tr (A* A A* A) 
+ A^tr (At A) + A 5 0t [At, A] 

-jL4> T ia 2 A^(t> + h.c. 



c. 

(5a) 

(5b) 



(5c) 



The covariant derivative of the Higgs triplet is given 
by 



D^A = d^A 



Hl-g 1 B fi A + ig 2 [W fi ,A] (6) 



and C is the charge conjugation matrix with respect to 
the Lorentz group. The counterterm parts of the La- 
grangian which are needed in the paper are given in 
App. [X] In addition, we consider an effective dimension 
5 (D5) operator which generates neutrino masses because 
it appears as soon as the Higgs triplet or a right-handed 
neutrino decouples: 



1 K ig {l i L UTi<j) (el C ta 2 ^j . (7) 



The most general neutrino mass matrix is given by the 
following formula 



(8) 



where k includes additional contributions to the dimen- 
sion 5 operator, like from gravitational effects [62j]. Thus 
the /3-function of the neutrino mass is given by the sum 
of the /3-functions for the contribution from the right- 
handed neutrinos and the contribution from the Higgs 
triplet. 



A 



A 



~ A 



FIG. 1: Feynman diagrams which are not considered in [6l(] 



The right-handed neutrinos and the Higgs triplet de- 
couple step by step at their respective mass scale and the 
effective theories have to be matched against each other. 
The decoupling of the right-handed neutrinos only con- 
tributes to the effective D5 operator. The decoupling 
of the Higgs triplet also gives a contribution to the SM 
model Higgs self-coupling because there is a coupling be- 
tween the SM Higgs doublet and the Higgs triplet given 
in Eq. (|5"c|) . Hence, the matching conditions of the right- 
handed neutrinos are 



.EFT 



2YjM~ 1 Y v 



and the decoupling of the Higgs triplet leads to 

Y A A 6 



.EFT 



A EFT = A + 2 



k — 2- 



Ml 

M 

Ml 



(9) 

(10a) 
(10b) 



In the following two sections, we present the calculation 
of the renormalization group equations. The calculation 
has been carefully double checked and compared to the 
results of Chao and Zhang [6l| . 



III. WAVE FUNCTION RENORMALIZATION 
AND COUNTERTERMS 

In order to obtain the wave function renormalization 
factor, all self-energy diagrams of the involved particles 
have to be calculated to one loop order. We use di- 
mensional regularization together with the MS scheme, 
because gauge invariance is generically preserved in 
this scheme. The differences to the formulae in [6l[ 
are underlined. Singly underlined terms differ in the 
prefactor and wavy underlined terms are not present 
in [6lj . Thus they probably stem from diagrams 
not taken into account in [6l| . For example, the 
wavy underlined terms in the counterterm of Ag are 
due to the Feynman diagrams shown in Fig[TJ Our 
Feynman rules and calculation can be downloaded from 



http : //www . mpi-hd . mpg . de/ ~mschmidt/rgeTriplet/ 

As we use a different definition for the couplings in the 
Lagrangian as in 61], the translation is presented in 
Tab. [T] The wave function renormalization factors are 
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our work 


Ya 


m 


Ma 


A 


Ai 


A 2 


A 4 


A 5 


A 6 


Chao and Zhang [61] 






M e 


A 




Ac 




-±A T 





TABLE I: Translation table for all relevant parameters 



given by 



SZ 



1 



A 



167r 2 e 



(3 - 6) 9\ + 4 (3 - 6) 52 " 2 tr y A ) 

(11a) 



5Z d 



1 



167r 2 e 
1 

167r 2 e 
1 

167r 2 e 



2T-^(3-ei).9?-^(3-6).9 2 2 



yjK + vty e + 3f2ya 



3 n 3 9 

YqCi.9i + 2^2.92 



(lib) 



(11c) 



where we have defined 

T = tr (yJy v + Y^Y e + 3YjY u + 3YjY a 



(12) 



The Yukawa coupling vertex £A£ can be renormalizcd 
multiplicatively 



1 



327r 2 e 



^(l-a)3? + (3-76) 52 2 



(13) 



the parameters in the Higgs potential, however, have to 
be renormalized additively 



Sm 2 



SMl 



6\ 
Mi 
<5A 2 

<5A 4 : 

SA 5 : 
SA 6 



1 



167r 2 e 
1 

167r 2 e 
1 

167r 2 e 
-8tr 
1 

167T 2 e 



3A - ^isi - ^ 2 <?2 j "' J - -i 0;: -^iv„) + <> \i - (i a, |J 

8A, + 2A 2 - ^i5? - -//i j i/1 + -I A ■///- + 2|A,,| 



6A 2 -1a(% 2 + 3 5: 



3 52 4 



Y^Y^Ye + YjY v YjY v + 3Y^Y U Y^Y U + 3YjY d YjY d 



12A? + 8A? 



12 \ () 70 

y63 2 ~ 8<& 2 J Ai + ^12 g 4 +18g 2 4 + - g 2 g 2 +14A 2 + 2A 2 + 4AiA 2 + 4 (A 2 + A 2 ) 



tr (YtY A YlY A 



1 



167r 2 e 
1 

167r 2 e 
1 

167r 2 e 



12 



144 



5 <i9l ~ 8. 9 ^ 2 J A2+I252 - —.91.92 + 3A^ + 12AiA 2 - 8A^+8tr 
3 ^ 2 n , 2 \, 27 4 



(yIy a yI 



(14a) 
(14b) 

(14c) 

(14d) 
(14e) 



-6gi + (8A1 + 2A 2 + 3A + 4A 4 ) A 4 + 8A 2 - 4 tr (yIy a Y^Y v 



\ti9l ~ y6ff 2 2 ) A 5 - y 5 2 g 2 2 + (2A X - 2A 2 + A + 8A 4 ) A 5 + 4tr (y a Y a Y}Y v ) 



1 



167r 2 e 



-1 h [yIyJmy v 



(14f) 
(14g) 

(14h) 



r 



IV. RENORMALIZATION GROUP EQUATIONS 

Using the vertex corrections and the wave function 
rcnormalization factors, we can deduce the /3-functions 
via the formula given in [r33| . As the wave function 
renormalization constant for the left-handed lepton dou- 
blets has an additional term with respect to the SM ex- 



tended by right-handed neutrinos, all vertices receive an 
additional contribution per left-handed lepton attached 
to the vertex which is given by 

\^ Y * (15) 
multiplied by the matrix characterizing the vertex from 
the left or from the right, respectively. In particular, the 
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/3-functions of the lepton Yukawa couplings [6J, [65|, \6t 



and the Yukawa coupling Y A become 



IS^Y V =Y V ^YJY„ - ^YjY e + ^Y A Y A 



Y fi 



T qi - -gi 

2CT 1 4 y2 

T - 4^1 - 4.92 



i67r 2 r A = 



1 ~y}y u + ~Y e ty e + ^y a y a y a + y a 



ii?y„ + iy e ty e + ^y a 



The renormalization of Ag is described by 



167T 2 A 6 



27 21 / \ 

* -T. 4A_4+_ 8As - — g\ - —g\ + 2T + tr (Y^ Y A J 



1CT 1 2 

and the anomalous dimensions of the Higgs triplet mass term is given by 



H>"'- w A =^ + 6^-4A 1 -A 2 -tr(Y2Y A 



(16a) 
(16b) 



Y A 



A fi - 2 tr 



-2A 4 m 2 - -|A 6 | 2 | M 



-2 
A 



The /3-function of the effective neutrino mass operator n function changes to 

-I T 



16tt 2 k 



K + K 



\yJy^ 3 -yX + IyIy a 



[2T - 3g 2 2 + A] k 



(16c) 



(17) 



(18) 



(19) 



r 



The RG equation for the type I contribution to the 
neutrino mass is only changed by the additional term to 
the /3-function of the neutrino Yukawa couplings due to 
the Higgs triplet. The remaining RG equations are pre- 
sented in App. [HI They either do not receive additional 
contributions or do not directly influence the neutrino 
mass matrix. The main difference in the RG equations 
compared to the results in [6ll | arise from contributions 
of the additional diagrams contributing to Aq. As it can 
been seen later, they have an impact on the evolution 
of neutrino masses, but the evolution of mixing angles 
and phases remains unchanged. In summary, the run- 
ning of the effective neutrino mass matrix m v above and 
between the seesaw scales is given by the running of the 
three different contributions to the neutrino mass matrix, 



C e Y}Y e + C„YjY u + C A YlY A 

(21) 



(i) 

where ml, stands for any of the three contributions to 
the neutrino mass matrix, respectively. The coefficients 
C e , v ,A and a are listed in Tab. [TTJ In the type-I+II see- 
saw scenario, large RG effects can be expected before the 
Higgs triplet is integrated out due to the different coef- 
ficients (C e , C Ul C A ) in analogy to the standard seesaw 
scenario where large RG corrections between the thresh- 
olds are induced by additional flavor-diagonal vertex cor- 
rections to the D5 operator. 



(i) v 

& = -y Y?M^Y V 
i& = f-A 6 M A 2 Y A . 



(20) 



The 1-loop /3-functions for m„ in the various effective 
theories can be summarized as 



16tt 



2 drriv 



(>) 



dt 



CeYlY e + C„Y}Y V + C A Y A Y A 



HIGGS TRIPLET IN THE MSSM 



In the MSSM, in addition to the Higgs triplet A ~ 
(3, 1), a second Higgs triplet A ~ (3,-1) with opposite 
hypercharge Y is needed to generate a D5 mass term 
for neutrinos. Furthermore, A ensures that the model is 
anomaly-free. Note, however, that only one Higgs triplet 
couples to the left-handed leptons. The additional terms 
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model 


M 
ml 


Ce 


c„ 


Ca 


flavor-trivial term a 


SM 




3 
2 


1 

2 


3 
2 


2T — 3g% + A 


SM 


2Y?M~ X Y V 


3 
2 


1 

2 


3 
2 


Z1 - TT)9l - 2^2 


SM 


-2A 6 M A 2 y A 


1 

2 


1 

2 


3 
2 


T - 2tr (Ia^a) - 3g| + A - 8Ai - 2A 2 - 4A 4 + 8A 5 - 
(4A 4 m 2 + |A 6 | 2 ) M A 2 - 4tr (y a Y?MY v ) Ag 1 


MSSM 


K 


1 


3 


3 
2 


2 tr (yjy„ + 3yty u ) + 3|A U | 2 - 2 (| 3 2 + 3 5 f) 


MSSM 


2YjM- 1 Y v 


1 


3 


3 
2 


2tr (yjy, + 3y„ f y„) + 3|a„| 2 - 2 (§<? 2 + 3<? 2 2 ) 


MSSM 


-2A 6 M A 2 y A 


1 


3 


3 
2 


2tr (yjy, + 3yjy) + 3|A„j 2 - 2 (f 5 2 + 3 5 f) 



TABLE II: Coefficients of the /3-functions of Eq. (|2f | , which govern the running of the effective neutrino mass matrix in 
minimal type II seesaw models. In the MSSM, the coefficients coincide due to the non-renormalization theorem [6?], [68| in 
supersymmetric theories. 





The decoupling of the Higgs triplet generates an effective 
dimension 4 term k eft in the superpotential, whereas 
the tree-level matching condition reads 



FIG. 2: Supergraphs which contribute to the wave function 
renormalization. 



.EFT 



' M A 



(23) 



in the superpotential are given by 



W A =M A Tr(AA) 



y/2 



V2 



V2 



(22) 



where D denotes the left-handed doublet and IhW denotes 
the Higgs doublets. We use the same notation as in [69| . 



The RG equations [8GJ can be obtained easily by using 
the supergraph technique as it is described in [69( . There 
are only two different types of supergraphs contributing 
to the wave function renormalization which are shown 
in Fig. [2] Here, we only show the RG equations which 
are relevant for the RG evolution of the neutrino mass 
matrix 



J 



167r 2 y A =Y A 



Y^Y e + Y,jY u + -YlY A 



YjY e +YjY v + ~Y A Y A 



Y A + Y A 



i|A d | 2 + itr(y^A)-^-7.g 2 2 



167T 2 A„ =A„ 



16tt 2 M a =M a 



2tr(3y u ty„ + YjY v ) + ^|A„| 2 - ^gf - 7g\ 
\ tr(YlY A ) + i|A tl | 2 + i|A,| 2 - 4 + 2g\ 



(24) 
(25) 

(26) 



As it can be seen in the /3-function of Y A for example, the sign of Y A Y A in P equals the one in the SM which leads 
to the same sign in the RG equations of the angles in the limit of a strong hierarchy (See Sec. IVI[) . The RG equation 
of the effective D5 operator k is 



I67T k =K 



y+y e + Y?n + _yt Ya 



Y^Y e + YjYu + -Y A Y A 



1 T 



K + K 



2tr(Y}Y v + 3Y^Y U ) + 3|A U | 2 - 2 ( -gf + 3 5 | 

(27) 
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There is also an additional contribution to the RG equations which are relevant for the type I contribution to neutrino 

masses 



16^y„ =Y V 



yty + 3yjy„ + -yIy a 



+ Y U 



tv(3Y^Y u + Y}Y V ) + ^|A U | 2 - (p 2 + 3.g 2 ) 



16n 2 M =M (2Y u Yj) T + (2Y v Yj) M 



(28) 
(29) 



All remaining RG equations are presented in App. [Cj 

The coefficients of the RG equations of the different 
contributions to the neutrino mass matrix are summa- 
rized in Tab. [TTJ Note, that the coefficients C e , C, 
Ca and a are the same for all three mass contributions 
rriu in the MSSM due to the non-renormalization theo- 
rem (sUsi. 



VI. RG EQUATIONS OF MIXING 
PARAMETERS IN THE TYPE II SEESAW CASE 

In order to understand the RG evolution of neutrino 
masses and leptonic mixing parameters in the presence 
of a Higgs triplet, we consider a type II model, where the 
neutrino mass is generated by a Higgs triplet only. The 
evolution of the mixing parameters in standard parame- 
terization can be described by the formulae of [52| with 
suitable replacements for P, F, a and a e : 



-P m^+m^P + am^ 



-YjY e ^YjY e + Y c !Y e F- 



a e YjY e 



(30a) 
(30b) 



Here, we can express P and F in terms of physical pa- 
rameters. 

P =C e diag(^,y2, y 2 ) + CaU * diag(y 2 ; y 2 ; y 2 )uT (31) 

F =D e diag(y 2 , y 2 y 2 ) + D A U* diag(y 2 , y 2 , y 2 )U T , 



where U is the MNS matrix and yi = ■ 
so-called standard-parameterization [7l| 



(32) 

We use the 



U = diag(e Ide , e 1 ^ , e ld ^)V diagte" 11 ^ 2 , e" IV2/2 , 1) 
V =# 23 (023 )r\Ri3 (e 13 )T s R 12 {fin) , 



(33) 



where Rij is the matrix of rotation in the i — j - plane 
and T.5 = diag(e i,5/2 , l,e~ iS/2 ). Note, that the Majorana 
phases drop out of the definition of P and F in flavor ba- 
sis. Hence the RG equations of the angles and the Dirac 
CP phase are independent of the Majorana phases, as it 
can be seen below. We derive the RG equations by using 
the technique described in the app endix of [52| which is 
based on earlier works [H, [zl, |73| . In the numerical ex- 
amples which are shown in the figures, we do not include 
any finite threshold corrections, since we are considering 
1 loop running and the finite threshold corrections are 



assumed to be of the order of 2 loop RG running. There- 
fore, the Higgs triplet is decoupled when its running mass 
equals the renormalization scale 



M A (/j dc 



(34) 



In all examples, we set Ma(Aqut) = 10 10 GeV. As we 
are only interested in showing the generic features of the 
RG evolution, we choose the Higgs self-couplings to be 
Ai,2,4,5 = 0.5 for simplicity, since they only indirectly 
influence the RG evolution of the angles and the flavor- 
dependent part of the RG equations of the masses. In a 
realistic model, the parameters Aj have to satisfy certain 
relations to produce the desired vevs. 

In the following, we present all formulae in the 
approximation y e <C y^ <C y T and 613 <C 
1. The exact formulae can be downloaded from 
http : //www.mpi-hd.mpg. de/~mschmidt/rgeTriplet. 



A. Running of the masses 

The main contributions to the RG equations of the 
masses 



16tt 2 — = Re a + 2C A ^r + 2C e y 2 sin 2 12 sin 2 (9 23 
mi v A 



0(0 



13 I 



(35a) 



16tt 



2 ™2 
7712 



Re a + 2C A - 
+ O(0is) 



2C e2 / 2 COS 2 012 Sin 2 023 



(35b) 



.2 m 3 
777 3 



16tt^ = Re a + 2C A ^r + 2C Wr cos2 ^3 + O(0 13 ) 

V A 



(35c) 



are the flavor-independent term Re a and the flavor- 

2 

dependent term 2 Ca ■ As the smallness of neutrino 
masses is usually explained by a small vev of the Higgs 
triplet va, the singular values yi = of the Yukawa 
coupling Ya can be of 0(1). This in turn leads to siz- 
able RG effects. Furthermore, the evolution of the mass 
squared difference is mainly given by 



A777 2 

16tt 2 — -£ w 2Rea + 4C A : '" J 
Ami 



(36) 



in the SM and MSSM with small tan/3. There can be a 
cancellation of the RG effect depending on the parame- 
ters A, in the Higgs potential and the sign of Ca, but 
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2 3 4 5 6 7 

(a)Evolution of neutrino masses 



9 10 11 12 13 14 15 16 
log |0 02/GeV) 



2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
log 10 0VOeV) 

(b) Evolution of the mass squared differences 



FIG. 3: As input values, we have chosen tribimaximal mixing at the GUT scale, mi = OeV, Am 2 tm = 2.5 10 3 eV 2 , Am^j = 
1.21(T 4 eV 2 , Ma = 10 10 GeV and A 6 = 2.510" 5 M A , corresponding to (A) = 0.15 eV. As we are only interested in showing 
the generic feature of the RG evolution, we choose the Higgs self-couplings to be Ai : 2,4,5 = 0.5 for simplicity, since they only 
indirectly influence the RG evolution of the angles and the flavor-dependent part of the RG equations of the masses. The 
shadowed area indicates the region where the Higgs triplet is present. It is integrated out at the border between the shadowed 
and the white area. 



generically the RG effect in the effective theory is large, 
as it can be seen in Fig. [3J This is just one possible ex- 
ample. The precise RG effect strongly depends on the 



parameters in the Higgs potential Aj. The charged lep- 
ton Yukawa couplings depend on the singular values of 
Yukawa coupling matrix Ya in a flavor non-diagonal way: 



J 



16tH 



to, 
, m 



■ Re a e + Da 



16tt 2 — tt =Rea e + D A 



16tt 



■ Re a e + Da 


















mi 




cos 2 


#12 " 




Wa 






V A 










'ml 






( ™l 


V 2 


sin 2 t 


'23 + 




- V A 






Wa 










'm\ 






( m 2 


k 


cos 2 ( 


?23 + 










Wa 



COS t>i2 



r 



0(6 13 ) 



7^ sin 1 



12 COS tf 2 3 



t-1 .9 

^ sm 2 ' 



12 sin 2 0- 



2:1 



0(9 13 ) 

D e y 2 T + 0(9 13 ) 



(37a) 
(37b) 
(37c) 



B. Running of the mixing angles 

Chao and Zhang [6l| have derived the formulae in the 
approximation |Y e | <C \Ya\ which captures the dominant 



effects as long as there is a strong hierarchy. Here, we cal- 
culate the renormalization group equations exactly (8lj 
and present the equations for the mixing angles in the 
approximation of vanishing y e , and 013: 



167T 2 0i2 = 



16tt 2 0i 3 = 



167T 2 02 3 = 



Am|i 2 (m 2 + mi) 

Da — 5 1- C e y T - 



v 



A 



sm V 23 



sin20 12 + 0(0 13 ) 



_e 2 (to 2 -TOi)to 3 

2 Vt (m 3 - mi)(m 3 - m 2 ) 



cos 5 sin 20i2 sin 20 23 + O (0 i3 ) 



; 'A 



6 1 • 9 

^ sin ( 



; 'A 



12 



■ COS 612 



2 m 2 - mim 2 + (m 2 - mi) m 3 cos20i2 
(m 3 - m 2 ) (m 3 - mi) 



(38a) 

(38b) 

sin 202 3 + (0i 3 ) 
(38c) 



where C e and Da are defined in Eqs. (|31l32p . The two neutrinos can be of the same order of magnitude and 
contributions to the running from charged leptons and 
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it strongly depends on the hierarchy of neutrino masses 
which of the two contributions is dominant. The con- 
tribution coming from the neutrino mass matrix (oc C e ) 
shows almost the same features as in the effective theory: 

• there is an enhancement factor which is propor- 
tional to t— \ , where mo denotes the mass scale of 
neutrinos; 

• the running strongly depends on tan/3 due to the 
charged lepton Yukawa couplings; 

• vanishing mixing is a fixed point. 

In contrast to the effective theory, however, there is no 
dependence on Majorana phases. This still holds for the 
exact equations. The RG evolution of the mixing angles 
is only influenced by the Dirac CP phase. On the other 
hand, the contribution from the charged leptons shows 
a completely different dependence on the Yukawa cou- 
plings, ft is basically proportional to the corresponding 
mass squared difference divided over the vev of the Higgs 
triplet squared. Hence, there is no large enhancement 
factor and no dependence on tan (3 in the SUSY case. 
Thus the overall size of the RG effect mainly depends on 
the vev of the Higgs triplet. 




• sin 26; 



(39) 



FIG. 4: Plot showing the evolution of the leptonic mixing 
angles in the SM. As input values, we have chosen tribi- 
maximal mixing at the GUT scale, mi = OeV, Amatm = 
2.510" 3 eV 2 , Am s 2 ol = 1.2f(T 4 eV 2 , M A = 10 10 GeV and 
A 6 = 2.5 10~ 5 Ma, corresponding to (A) = 0.15 eV. As we are 
only interested in showing the generic feature of the RG evo- 
lution, we choose the Higgs self-couplings to be Ai,2,4,s = 0.5 
for simplicity (In a realistic model, the parameters A, have to 
satisfy certain relations to produce the desired vevs.), since 
they only indirectly influence the RG evolution of the an- 
gles and the flavor-dependent part of the RG equations of the 
masses. In the MSSM with small tan/3, the plot looks similar 
due to Ca being the same in both theories. The shadowed 
area indicates the region where the Higgs triplet is present. 
It is integrated out at the border between the shadowed and 
the white area. 



This gives a good estimate for the running in the strongly 
hierarchical case. The sign of the RG effect is determined 
by the sign of the mass squared difference and the factor 
Da in front of the factor Y^Ya in P. As Da is positive 
in the SM and MSSM, 823 is evolving to larger values 
coming from the high renormalization scale for a normal 
hierarchy. Furthermore, the /3-function is approximately 
proportional to sin 26ij which implies that a vanishing 
angle remains small. Taking into account these generic 
features, the RG effect from the charged leptons is largest 
on #23 due to the combination of a large mass squared 
difference and a large angle. As it can be seen from the 
equations, zero mixing is a fixed point. This is also ob- 
vious from the RG equation in matrix form: in this con- 
figuration, P and F will be diagonal, if Y e and Y v are 
diagonal. In FigfJl we have plotted the evolution of mix- 
ing angles in the SM for a strongly hierarchical spectrum 
in order to suppress the effect coming from the effective 
D5 operator. The gross features of the running can be 
immediately seen: the only sizable effect is on 623 due to 
the large angle and mass squared difference. As it can be 
seen from the plot, the RG effect can be estimated by a 
leading log approximation to 



is given by 



D A 



^ ( ^ - ^ cos 2 9 12 - ^ sin 2 0,2 ) sin 29 13 (41) 



rn. 



nr, 



Let us comment on the configuration 613 — 7x13 = 0, 
which is stable under the RG in the effective theory. Van- 
ishing mass eigenvalues remain zero, as it can be seen 
from Eq. (|35cp . but 613 receives corrections 



16tt 2 0i 



' 2 «A (Vr - vl) (vl - vt) 

x sin 20 23 +0(013, y 3 ) 



cos S sin 26*12 



(42) 



Thus #13 = 7713 = is not stable under the RG. However, 
the effect is negligible, because (jf^j ^^a 1 ' is very small 
and m 3 — is stable. 



A0, 



Da Ar 



— sin 26L In 

3 Ma 



(40) 



The contribution to 613 coming from the charged leptons 
vanishes in our approximation. For non-vanishing #13, it 



C. Running of the phases 



The RG evolution of the phases is rather small: 
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16tH(5 



167T Ifii 



167T 02 



C (m2 mi)m3 sin J sin 20 12 sin 20 23 era 1 +0(8 13 ) 



2 (m 3 - mi)(m 3 - m 2 )' 



(43a) 



2 {m\ + m|)m 2 sin 2 6*23 — ((m? + m 2 ) sin 2 612 — mim2(cos 20i2 — cos 2# 2 3))to3 



(ni3 — mi)(m3 — TO 2 )(m 2 — mi) 



cot( 



D A - 



'■21 



2C e y 2 T - 



sin 2^i2 



cot 23 sm5 13 + O{9 2 13 ) 



(43b) 



„ Am 2 , 
■Da— 2^8^2012 



-(m? + m\)mj, cos 2 #12 + mi((m 2 + m|) sin 2 023 + m 2 m 3 (cos 20 i2 + cos 2023)) 
(m 3 - mi)(m 3 - m 2 )(m 2 - mi) 

cot6>23sm<5 0i3 + 0(0 2 3 ) , 

I 



tan 0i2 



(43c) 



because the leading order of the Majorana phases is of 
order 813. Only the Dirac CP phase S involves a term 
which is inversely proportional to #13. Thus, there is a 
sizable effect for small #13. For vanishing 0i3, 6 has to 
vanish (for realistic values of O12 and 023) in order to 
ensure analyticity of S(t). The RG equation of the Dirac 
CP phase S does also not depend on the Majorana phases 



VII. RG EVOLUTION IN THE FULL TYPE II 
SEESAW CASE 

In the full type II case, it is not possible to ex- 
press the RG equations in terms of mixing parame- 
ters. Therefore one has to resort to numerical calcula- 
tions. For this purpose, we have extended the Mathe- 
matica package REAP, which is available on the web page 



http : //www . ph . turn . d e/ Tge[ to include a left-handed 
triplet. 

To illustrate the largeness of RG effects in the full type 

11 seesaw scenario, we show an example, where bimaxi- 
mal mixing at high energy evolves to the LMA solution 
at low energy. In previous works [3(| HH, [H, H3] , this evo- 
lution was due to an inverted hierarchy in the neutrino 
Yukawa couplings Y v or large imaginary off-diagonal en- 
tries. Here, the relevant matrix YJY V is real and has 
a normal hierarchy. In addition, the singular values of 
the Yukawa coupling matrix Y"a are small (0(1O -5 )). In 
spite of the small couplings, there is a sizable effect on 

012 which can be seen in Fig. [5l It is due to the different 
RG equations of the contributions to the neutrino mass 
matrix. 

In our example, we have chosen Aq to be relatively 
large A 6 = O(10 9 ) GeV, because it receives corrections 
of the order of M3 (Y„) 2 3 (1a)33< The evolution of the 
mixing angles 0i2 and 023 is highly non-linear above the 
threshold of the Higgs triplet. Hence, a leading log ap- 
proximation is not possible. In the MSSM, the equations 
for the mixing angles presented in [52j are valid at each 



renormalization scale /i. Hence, 0i2 is increasing, as long 
as there there are no imaginary off-diagonal entries and 
there is a normal hierarchy in the neutrino Yukawa cou- 
plings. 



VIII. CONCLUSIONS 

We calculated the RG equations in the type II seesaw 
case and found differences to the ones calculated by Chao 
and Zhang [6l[ in the RG equations of the parameters of 
the Higgs potential. 

In the SM, the matrix P describing the off-diagonal 
contributions to the RG evolution of the neutrino mass 
matrix is different for the contribution coming from the 
Higgs triplet compared to the one for the effective D5 
operator. Hence, there can be large RG effects due to 
the different RG equations of the different contributions 
to the neutrino mass matrix. 

Furthermore, we derived the exact RG equations in 
terms of the mixing parameters. The equations have 
a different structure compared to the ones in the stan- 
dard seesaw case as well as in the effective theory case. 
One feature is that the RG equations of the mixing an- 
gles are independent of the Majorana phases. There- 
fore, there is no damping of the RG effect due to Majo- 
rana phases. The main difference to the running in the 
standard seesaw scenario is the proportionality of the (5- 
function of the mixing angles to the mass squared dif- 
ference in contrast to the inverse proportionality in the 
case of a hierarchical spectrum. Hence, there is no en- 
hancement factor and the RG effect is small as long as 
Ya is small. Furthermore, the RG effect of the mixing 
angles 0^ is proportional to sin 20^ which leads together 
with the proportionality to the mass squared difference 
to a large effect on 023 compared to the effect on 0i2 
due to Am 2 ol <C Am 2 tm and compared to 0i3 due to 
sin 20i3 < sin 20 23 . 

The RG equations in the full case can only be stud- 
ied numerically. The interplay of the contributions from 




2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
log 10 0i/OeV) 

(a)Angles 




2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
log 10 MGeV) 

(b)Masses 



FIG. 5: In the full type II seesaw case, there is a compli- 
cated interplay between the two contributions to the neu- 
trino mass matrix. Here, we just plot an example for the 
following initial values at the GUT scale: M A = 10 10 GeV, 
A 6 = 4.56 10 9 GeV, mi = 0.02 eV, Am s 2 d = 1.5 10" 4 eV 2 , 
AmL = 5.5 l(T 3 eV 2 , 9 12 = 823 = f , 013 = 0, 5 = ^ = 
V2 =0,Y U = 0.37 diag (l0~ 2 , 10~\ l), where Y A is chosen 
diagonal Y A = diag (1.3 10" 5 , 1.5 10" 5 , 5.1 10~ 5 ) and M is 
chosen appropriately to produce bimaximal mixing. The dif- 
ferently shaded areas indicate the different energy ranges of 
the various effective field theories. At each border, a par- 
ticle, either a right-handed neutrino or the Higgs triplet, is 
integrated out. 
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APPENDIX A: COUNTERTERM LAGRANGIAN 

The relevant wave function renormalization factors are 
defined in the usual way by 

(£ L ) B =ZfJ L (Al) 
A B =Z A A (A2) 
cb B =z\$ . (A3) 

The Yukawa couplings are renormalized multiplicatively 

{{ZJ L y (Y A ) B zf L Zi) fg =/x* (YAZ Y J fg (A4) 
(zl R {Y V ) B Zl Z l) fg =^ ( Y » Z »)fg ■ ( A5 ) 

The parameters in the Higgs potential have to be renor- 
malized additively 



Z,pm B 


—m 2 + 8m 2 


(A6a) 


Za(M 2 a ) b 


=M A + SMl 


(A6b) 




=^Z X X 


(A6c) 


Zl (A 1>2 ) fl 


=M £ (Ai l2 + *Ai, 2 ) 


(A6d) 


Z A Z,p (A 4i5 ) B 


=tf (A 4:5 + £A 4 ,s) 


(A6e) 


z\y {zif {u) B z\ 


=H* (A 6 + <5A 6 ) . 


(A6f) 



right-handed neutrinos and the Higgs triplet can lead to 
large RG effects even in the SM. Hence, it is necessary to 
consider RG effects in model building to make predictions 
which can be compared to the experimental data. 



The insertion of the definitions of the renormalized quan- 
tities into the bare Lagrangian yields the counterterm 
part of the Lagrangian which is needed to cancel the di- 
vergences. Thus the counterterm Lagrangian 



Acknowledgements ^ = ^ + Ca + Cy(A) (A7) 

M.S. wants to thank Werner Rodejohann for useful dis- 
cussions and proofreading the manuscript, as well as He is given by 



(SZ VR ) 



•■"'at ''' 2''/.''' 



{5Z M M).v f R - {Y v 8Z Yu ) qf V^ 



>gf "R 



gf 1 



(SZ VR ) gf P K 



h.c. 



(A8a) 
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C A =SZ A tr(D^(D^A) + \iti(D^ (sZ gi9l B»A + 5Z g2 g 2 ^ [W^,A]) + h. 



C V (A) =SM A tr A*A + ^ (tr A f A) 2 



SA 2 
2 



(tr A^) 2 - tr (A^AA^A) 



6A 6 



SA^UtrA^A + SA^ [A\ A] + ^4> T (ia 2 ) A\ 



(A8b) 



(A8c) 



APPENDIX B: RG EQUATIONS IN THE SM 

The remaining RG equations of Yukawa coupling matrices 



16n 2 Y d 


=Y d 


2 d d 


-yty 

2 " u 


+ y 


16tt 2 Y u 


=Y U 


-Y^Y - 

2 u u 


2 d d 


+ y« 



where T = Tr 



y e ty + yjy + 3yjy d + 3y]y 



and the right-handed neutrino mass matrix 

T 



i6tt 2 m = (y,yj) M + M (Y u Yjy 



are listed for completeness. They are taken from |52j ]. The anomalous dimensions of the Higgs doublet mass 
by 
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9 



A/ 



,|A 6 | 2 



16*%» = Yq9 1 + ~/ 2 ~ ~ 2 X - T + - 2 tr (yjM »y„ ) - 3A 4 -A - 3 ^ 
and the renormalization of the remaining couplings in the Higgs potential are described by 

2 

i / 3 2 i 2 ^ 
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I(»,t-A =(iA- :!,\ ( 3g 2 + |g 2 J + 3<?f + | Q<7? • .</ I f 4AT 



16tt 2 Ai 



16tt 2 Ao = 



16tt 2 A 4 = 



167T 2 A 5 



8 tr (y^yyjye + yjyyjy + 3yjy„yjy„ + syJy^y^ + 12A 2 + sa 2 

f 5 2 A X - 24 5a a A! + ^ 5l + 18 52 4 + ?3 2 3 2 2 +14A 2 + 4AiA 2 + 2A 2 + 4A 2 + 4A 2 
5 25 5 

4tr (y A y A ) A x -8tr (y A y A y A yA 



y ff 2 A 2 - 24 ff 2 A 2 +12gJ - ^g 2 g 2 2 + 3A 2 + 12AiA 2 - 8A 2 + 4tr (V A y A ) A , 



tr (yIy a yIy a 



-<?? A4 - y3 2 2 A 4 + ^gf+6 g: 



8A1 + 2A 2 + 3A + 4A 4 + 2T + 2 tr y A y A 



A 4 



8A 2 -4tr [YlY A Y}Y v 
L 2 A 5 -^ ff 2 A 5 18 ~ 2 - 2 



2Ai - 2A 2 + A + 8A 4 +2T + 2 tr 



(yIy a ) 



A s 



4tr 



(y A y A yjy, 



There is also a contribution to the gauge coupling renormalization due to the Higgs triplet 



41 



1 



.1 1 47 
16^ 2 /3 S1 =77-5? + - • 3 • 2-g\ = -g* 
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G 



16tt 2 /3, 



<J2 



19 , 1 

On H 2-2(7, 

6 y2 6 y2 



5" 1 10 

3 



5 3 
2 52 



(Bl) 
(B2) 

(B3) 
is given 

(B4) 



(B5a) 
(B5b) 

(B5c) 
(B5d) 

(B5e) 

(B6) 
(B7) 
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which has been also derived in [6l[. The RG equation of the strong coupling constant remains unchanged 

16n 2 f3 g3 = -7gl . 



(B8) 



APPENDIX C: RG EQUATIONS IN THE MSSM 

The RG equations which do not directly influence the renormalization of the operator generating neutrino masses 
are given by: 



16tt z Y k =Y k 



l6ir 2 Y u =Y U 



3Y*Y e + Y}Y V + -Y^Ya 



16ir 2 Y d =Y d \3YjY d + Y^Y U 



Yd 



Y e tr(3Y^Y d + Y}Y e ) + ^\A d \ 2 - ^gf - 3g: 
tr(3rJy M + Yfr,) + 3 -\A u \ 2 ^gf - 3g 2 2 - fg 2 3 



tr{wlY d + r e ty e ) + ^|A d | 2 - L 9 i _ 3g 2 _ 



1G 



167r 2 A d =A d 



2tr(3Y^Y d + Y}Y e ) + 7 -\A d \ 2 + 1 tr(YlY A ) - ^g 2 - 7g 2 



(CI) 
(C2) 
(C3) 
(C4) 



There are also contributions to the evolution of the gauge couplings of electroweak interactions: 

16tt 2 /3 Si =^gf + 2 • 3^gf = ^gf 
16vr 2 /3 92 -lg 2 3 + 2-2 32 3 = 5 ff2 3 . 
The RG equation of the strong coupling constant remains unchanged 



(C5) 
(C6) 



16n 2 p ga = -3g 3 3 



(C7) 
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